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ABSTRACT 

In the past few years, several disks with inner holes that are relatively empty of small dust grains 
have been detected and are known as transitional disks. Recently, Spitzer has identified a new class 
of "pre-transitional disks" with gaps based on near-infrared photometry and mid-infrared spectra; 
these objects have an optically thick inner disk separated from an optically thick outer disk by an 
optically thin disk gap. A near-infrared spectrum provided the first confirmation of a gap in the 
pre-transitional disk of LkCa 15 by verifying that the near-infrared excess emission in this object 
was due to an optically thick inner disk. Here we investigate the difference between the nature of 
the inner regions of transitional and pre-transitional disks using the same veiling-based technique to 
extract the near-infrared excess emission above the stellar photosphere. However, in this work we use 
detailed disk models to fit the excess continua as opposed to the simple blackbody fits used previously. 
We show that the near-infrared excess emission of the previously identified pre-transitional disks of 
LkCa 15 and UX Tau A in the Taurus cloud as well as the newly identified pre-transitional disk of 
Rox 44 in Ophiuchus can be fit with an inner disk wall located at the dust destruction radius. We also 
present detailed modeling of the broad-band spectral energy distributions of these objects, taking into 
account the effect of shadowing by the inner disk on the outer disk, but considering the finite size of the 
star, unlike other recent treatments. The near-infrared excess continua of these three pre-transitional 
disks, which can be explained by optically thick inner disks, are significantly different from that of the 
transitional disks of CM Aur, whose near-infrared excess continuum can be reproduced by emission 
from sub- micron-sized optically thin dust, and DM Tau, whose near-infrared spectrum is consistent 
with a disk hole that is relatively free of small dust. The structure of pre-transitional disks may be a 
sign of young planets forming in these disks and future studies of pre-transitional disks will provide 
constraints to aid in theoretical modeling of planet formation. 

Subject headings: accretion disks, stars: circumstellar matter, planetary systems: protoplanetary disks, 
stars: formation, stars: pre-main sequence 



1. INTRODUCTION 

Several disks which have nearly photospheric near- 
infrared emission but substantial excesses above the stel- 
lar photosphere at wavelengths beyond ~20 /zm have 
been observed and are referred to as "transitional disks" 
(|Strom et al.1 IT989). Usin g data from the S pitzer In- 
frared Spectrograph (IRS: iHouck et all 12004V detailed 
modeling has demonstrated that this flux deficit at near- 
infrared wavelengths relative to full disks can be ex- 
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plained by optically thick disks with inner holes of < 
40 AU. In most cases these inner holes are not com- 
pletely devoid of material (e.g. GM Aur, TW Hya, 
CS Cha, CVSO 224); a minute amount of micron- or 
sub-micron-sized optically thin dust exists within the 
hole, producing a small infrared excess over the photo- 
spheric flux, still well below the med ian excess of Class 
II ob j ects, as well as silicate emissio n (|Calvet et al.ll2002l 
[20051 lEspaillat et all |2007pI l2008bl ). Gas has also been 
detected within the inner holes of tra nsitional dis ks (e.g . 
iNaiita et all [2003: Ber gin et al]l2004t iSalvk et al.ir2 007). 

Recently, the Spitzer Space Telescope (Werner et al. 
I2004f ) identified a new class of disks called "pre- 
transitional disks" ar ound LkCa 15 and UX Tau A 
IjEspaillat et al. 2007a|). These disks have deficits of mid- 
infrared flux (5-20 fim) and substantial excesses at longer 
wavelengths, as is seen in the transitional disks. How- 
ever, in contrast to the small or absent near-infrared 
(2-5 /mi) excesses exhibited by transitional disks, pre- 
transitional disks have significant near-infrared excesses 
relative to their stellar photospheres, similar to the me- 
dian sp_e£tral_encrgy distribution (SED) of disks in Tau- 
rus (jD'Alessio et al . 1999). The distinctive shapes of 
these SEDs indicate that pre-transitional disks have an 
inner disk separated from an outer disk and that we may 
be seeing the development of gaps within protoplanetary 
disks. 

While the truncation of LkCa 15's outer disk has 
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been imaged in the millimeter (Pi etu et al.l I2006T) , 
lEspaillat et al.l (|2007al ) showed that the substantial near- 
infrared excess of LkCa 15 could be explained by either 
optically thick material or by ^10~ n M© of optically 
thin dust mixed with the gas in the i nner di sk. In order 
to resolve this issue, lEspaillat et al.l (|2008a[ ) obtained a 
medium resolution near-infrared spectrum spanning the 
wavelength range 2-5 /im. This near-infrared spectrum 
had absorption lines that were weaker relative to the 
spectrum of a standard star of the same spectral type. 
This phenomenon, known as "veiling" (Hartigan et al.l 
1989), is also observed in similar spectra of full disks and 
is due to emissi on from dust located at the dust sub- 
limatio n radius (|Muzerolle et al.l [20031) . lEspailla t et al.l 
(|2008af) measured a veiling factor (r^-) of 0.3±0.2 for 
LkCa 15 at ^2.2 /im and fit the near-infrared excess from 
2-5 /im with a single-temperature blackbody of 1600 K. 
This behavior can be explained by an optically thick wall 
locate d at the dust sublima tion radius, as is seen in full 
disks (Muzcrol le et al.| [2003) . These data confirmed that 
LkCa 15 has an inner optically thick disk, making this 
observation the first independent verification of a gap in 
a protoplanetary disk. 

Here we expand our sample to include the pre- 
transitional d isks of UX Tau A ( jEspaillat et al.ll2007al) 
and Rox 44 ()Furlan et al.l I2009D , and the transitio nal 
disks of GM Aur and DM Tau (jCalvet et al.l 120051) in 
order to explore the structure of the inner regions of 
pre-transitional and transitional disks. To do this we 
use veiling measurements in the K-band to extract the 
near-infrared (-IR) excess emission of t hese objects and 
then fit this emission w ith disk models (D'Alcssi o et al.l 
120051: iCalvet et al.ll2002t ). Given the sensitivity of veiling 
measurements on the adopted spectral type, we redeter- 
mined spectral types for our targets using optical spectra 
(§ 14. ip . We also tested our veiling measurement methods 
with the spectrum of the diskless, weak-line T Tauri star 
(WTTS) LkCa 14 rS HTB . 

We find that the near-infrared spectra of the pre- 
transitional disks of LkCa 15, UX Tau A, and Rox 44 
are well-explained by the wall of an optically thick inner 
disk (§ |L2~T1 S lLiOl § I4T2T4]) . In contrast, our data shows 
that the inner hole of the transitional disk of GM Aur 
contains a small amount of optically thin sub-micron- 
sized dust while DM Tau's hole is relatively free of small 
dust (S 14.2.51 fc § I4.2.6[) . Our results are consistent with 
veiling and interferometric measurements found in the 
literature. We also perform detailed model fits to the 
broad-band SEDs of our pre-transitional disk sample, 
and explore the effect of shadowing of the outer disk 
by the inner disk (§ 14.31) . taking into account that the 
star is not a point source, as has been adopted in other 
studies (|Espaillat et al.ll2007al : iMulders et al.ll2010[ ). The 
structure of pre-transitional disks suggests that of the 
disk cleari ng mechanisms proposed to date, planet forma- 
tion (e.g. iGoldreich fc Tremainel fl980l iRice et aT1l200l 
iVarniere et al.l 120061) is most likely a dominant factor in 
clearing these disks. 

2. OBSERVATIONS & DATA REDUCTION 

Near-infrared spectra of LkCa 14, LkCa 15, UX Tau A, 
Rox 44, GM Aur, and DM Tau were obtained at the 
NASA Infrared Telescope Facility (IRTF) facility using 
SpeX (jRavner et al.l 120031 ) . For each of our targets, we 



used the Long XD2.1 grating, covering 2.1 to 5.0 /Ltm, 
with a slit of 0'.'5 x 15" for a resolution (X/SX) of 1500. To 
trace near-IR emission at shorter wavelengths, if present, 
we obtained additional spectra from ~0.8-2.5 /im for 
LkCa 15, UX Tau A, Rox 44, and GM Aur using the 
low-resolution prism with a slit of 0'.'8xl5" (A/<5A^250) 
with an exposure time of 40 seconds. The dates of the 
observations, the LXD exposure times, and the signal- 
to-noise ratios for each target are given in Table [TJ We 
n ote that the prism dat a for Rox 44 were first presented 
in McClureet all (|2010t l. 

We extracted the data with Spextool (jCushing et al.l 
I2004T) and co rrected for telluric absorption using the xtell- 
cor routine (jVacca et al.l l2003h . Our telluric standards 
for LkCa 14, LkCa 15, DM Tau, GM Aur, UX Tau A, and 
Rox 44 were HD25175, HD27777, HD27761, HD27777, 
HD27777, and HD146606 respectively. Bad pixels and 
regions of high telluric noise (e.g. 2.5-2.9 ^m and 4.2- 
4.6 /im) were manually removed from the spectra. 

To measure the spectral types for objects in our sam- 
ple, we used archival low-dispersion optical spectra ob- 
tained at the 1.5 m telescope of the Whipple Obser- 
vatory with the F ast Spectrograph for the Tillinghast 
Telescope (FAST; iFabricant et ~atlll998l ) equipped with 
the Loral 512x2688 CCD. The spectrograph was set up 
in the standard configuration used for FAST COMBO 
projects: a 300 groove mm -1 grating and a 3" wide slit. 
This combination offers spectral coverage over ^3600- 
7500 A, with a resolution of ^6 A. All spectra were 
reduced at the Harvard- Smithsonian Center for Astro- 
physics using software developed specifically for FAST 
COMBO observations and were wavelength-calibrated 
and combined using standard IRAF routines. Data for 
UX Tau A and LkCa 14 were obtained in 1995 and 1996 
as part of Program 38 (PI: Briceho) and are publicly 
available in the FAST database03 Spectra for LkCa 15, 
DM Tau, and GM Aur were obtained in 1995 and 1996 as 
part of Program 30 (PI: Kenyon) and were first used in 
iKenvon et al.l (|1998l ) to measure the equivalent widths of 
H a , [N II], He I, and [S II] lines. We note that no FAST 
spectra were available for Rox 44. 

3. DISK MODEL 

In the next section, we fit the near-IR excess emission 
and the broad-band SEDs of some of the objects in our 
sample with disk models. Below we describe the models 
used to reproduce the observed emission. 

In the pre-transitional disks, there is an inner wall lo- 
cated at the dust sublimation radius which dominates 
the near -IR (2-5 urn) emission . Here we use the wall 
model of lD'Alessio et al.l (|2005| ). The wall is assumed to 
be verticaO with evenly distributed dust. The wall's op- 
tical depth increases radially and is optically thin closest 
to the star. The stellar radiation impinges directly onto 
the wall, with an angle of 0° to the normal o f the wall's 
surface. Following ICalvet et al.l (J199U 11992t h the radial 
distribution of temperature for the wall atmosphere is 
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10 http://tdc-www.harvard.edu/cgi-bin/arc/fsearch 

11 We leave it to future work to explore the curvature of the 
inner wall (Nagel et al. in prep). 
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where F Q = (X* + L acc )/-iTTR 2 waa , <r R is the Stefan- 
Boltzmann constant, T4 is the total mean optical depth 
in the disk, and n s and k^ are the mean opacities to the 
incident and local radiation, respectively. L* is the stel- 
lar luminosity and L acc (~ GM*M / R*) is the luminosity 
of the accretion shock onto the stellar surface, which also 
heats the inner wall. To derive the location of the wall, 
one sets t<j = and 
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T W aii is the temperature at the surface of the opti- 
cally thin wall atmosphere (i.e. Td( T d—0))- When fit- 
ting the wall located at the dust sublimation radius, 
the relevant dust grains located at these high temper- 
atures (1000-2000 K) are co mposed of silicates. Fol- 
lowing D'Ales sTo et al.l (|2005l) . here we adopt silicates 
with a dus t-to-gas mass rat i o (Cif ) of 0.0034 using opac- 
ities from iDorschner et al.1 (|1995l ). The grain size dis- 
tribution used here (and throughout the disk) follows 
the form Q~ 3 ' 5 whe re a varies between a m i n and a max 
(jMathis et al.lfl977h . The temperature of the inner wall 
(1 wail) > the height of the wall {z l wall ), and the maximum 
grain size are adjusted to fit the SED. 

In the transitional disks, there is an inner hole that 
sometimes contains a small amount of optically thin dust 
which contributes to the SED between ~2-10 /im. We 
calculate the emission from this optically thin dust region 
by integrating the flux from optically thin annuli where 
all of t he dust grain s are h eated by stellar radiation fol- 
lowing iCalvet et al.l (|2002h . Some pre-transitional disks 
also contain optically thin dust within their disk gaps. 

In both pre-transitional and transitional disks, there 
is an outer wall located where the outer disk is trun- 
cated and this wall dominates the SED emission from 
~20-30 /im. The temperature structure of this wall is 
calculated in the same manner as the structure of the 
inner wall above. For the wall of the outer disk, we 
ad d organics and t roilite to the dust mixture follow- 
ing |DiAlessioIeFaI] ([2005| ) with ( org = 0.001 and C, tro i = 
0.000768 and sublimation temperatures oiT org = 425 K 
and Tt r oi = 680 K. We include water ice as well where 
the local temperature is below the sublimation temper- 
ature of this component (Ti ce — 110 K). We use Q ce = 
0.00056, which is 1 % of the a bundance proposed by 
iPollacket all (119941) . sinc e ID'Alessio efaD (|2006f ) found 



that [Pollack et al. (199J)'s ice abundance produces fea- 
tures that are not observed in typical disk SEDs. Opac- 
ities for organics, troilit e, and water ice a re ad opted 
from iPpllack et al.l (J1994D . iBegemann et al.l ()1994D . and 
iWarren (1984). In this paper we do not study the de- 
tailed composition of the dust, but simply illustrate that 
a typical dust composition can reasonably explain the ob- 
served SED. We adopt an ISM-sized grai n distribution 
(i.e. a m i n — 0.005 /im and a rna:r =0.25/im; iMathis et al.l 
|1977|) and varied the outer wall temperature (T^ a;/ ) and 
height (z° aH ) to achieve the best- fit to the SED. 

In the pre-transitional disks, the inner optically thick 
disk will cast a shadow on the outer disk. Assuming that 
the central star is a point source, this shadowing of the 
outer disk would be substantial. On the contrary, the 
star appears as a finite source to the outer disk and so 
the inner disk will cast a large penumbra and a small 



umbra on the outer disk. If a point on the outer wall 
is in the umbra, it does not see the star. If it is in the 
penumbra, it will see some of the star, and will still be 
illuminated to some extent. If the point on the wall is 
completely out of the shadow (both umbra and penum- 
bra), it is fully illuminated by the star. In calculating the 
outer wall's emission we take into account shadowing by 
a finite source star (see Appendix). Since we assume the 
wall is vertical and the luminosity along the the outer 
wall changes due to the geometry of the shadow, there 
will be a range of temperatures along the surface of the 
wall, consistent with the varying degree of illumination. 
We calculate the resulting SED with this range of tem- 
peratures and the radius of the wall is determined by the 
temperature in the fully illuminated part. 

Both pre-transitional and transitional disks have a con- 
tribution to their SEDs beyond ~40 /im from the outer 
disk and we use the i rradiated accretion disk models of 
ID'Alessio et al.l (|2006f ) to model the emission of the disk 
behind the wall. To simulate grain growth and settling, 
the disk is composed of two dust grain size distribu- 
tions. In the upper disk layers, grains are ISM-sized and 
in the disk midplane the maximum grain size is 1 mm 
(jD'Alessio et al.l[2006| ). The viscosity parameter (a) and 
settling parameter (e; i.e. the dust-to-gas mass ratio in 
the upper disk layers relative to the standard dust-to-gas 
mass ratio) are varied to achieve the best-fit to the SED. 
We adopt an outer disk radius of 300 AU for all of our 
disks. 

4. ANALYSIS 

Classical T Tauri stars tend to have observed pho- 
tosphcric absorption lines which are weaker than those 
of non-accreting stars. This "veiling" is due to an ex- 
cess continuum which adds to the star's emission and 
"fills in" stellar absorption lines, making them appear 
weaker when compared to a stand ard star of the same 
spectral type (|Hartigan et al.l!T989l ). In the optical and 
shorter wavelengths, this is primarily due to the hot con- 
tinu um which arises from the accretion shocks onto the 
star (jCalvet fc Gullbringlll998l ). In the near- infrared, for 
stars which are not strong accretors, veiling of photo- 
spheric absorption line s is mainly due to dust in the disk 
(|Muzerolle et all 12001 . 

Since veiling measurements are dependent on the spec- 
tral types that are adopted, we first re-derive spectral 
types for most of our sample using FAST optical spectra 
in § 14.11 Then, in § !4.2l we measure near infrared veilings 
for our sample, extract the excess emission, and fit it with 
disk wall models. We then model the broad-band SEDs 
of the pre-transitional disks of LkCa 15, UX Tau A, and 
Rox 44 with the best-fit inner wall model, including the 
shadowing of this inner wall on the outer disk. 

4.1. Spectral Type Classification 

Spectral types were derived with the SPTCLASS 
toolQ an IR AF/IDL code based on the methods de- 
scribed by iHernandez et al.l ()2004h . The original code 
was optimized to calculate the spectral types of early 
type stars (B, A, & F) based on 33 spectral features that 
are sensitive to changes in T B g. SPTCLASS was later 

12 http://www.astro.lsa.umich.edu/~hernandj/SPTclass/ 
sptclass.html 
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revised to incorporate 11 spectral indices optimized for 
solar type stars (from middle F to early K) and 16 spec- 
tral indices optimized for low mass stars (from early K 
to M5). Each spectral index is obtained by measuring 
the decrease in flux from what would be expected at the 
line center by interpolating between two adjacent bands 
on each side of the line. Each index is calibrated using 
spectroscopic standards observed with FAST. 

Most of the spectral classifications derived in this work 
differ somewhat from those reported in the literature. 
Our classification of DM Tau (M1.5±1.0) is in agreement 
with previous measurement s of MO. 5 (jCohen fe Kuhil 
I1979D and M2 (|Herbidli"977l ). However, we measure a 
spectr al type for LkCa 14 of K5.5±1.0 while iHerbig et all 
(|1986j) derived a spectral type of MO. For GM Aur, we 
derive a spectral type of K5.5±1.0, but in the litera- 
ture we find spectral cla ssifications of K3 (Hcrbig 119771 ) 
and K7-M0 IjCohen fe K uhi 1971) for this object. We 
derive a spectral type of K3±1.0 for LkCa 15 which dif- 
fers from__the_K5 spect ral type for this object reported 
by IHerbig et all (119861) W e cla ssify UX Tau A as a 
G8.0±2.0 star; IHerbig] (QIZZD and lCohen fc Kuhil (11971) 
both reported a spectr al type of K2, Rvdgrc n et al.l 
(|1976t l reported G5, and lHartigan et al.l (|1994h find K5- 
K2. Since no FAST spectra are available for Rox 44, we 
adopt a spectral type of K3 from lBouvier fe Appenzellerl 
(J19921) for this object. 

The spectral types derived here (Table [2]) are based 
on a greater number of spectral indices than used in the 
works cited above. Since commonly used spectral indices 
such as He I (4922 A, 5016 A) and Na I (5890 A, 5896 A) 
can be affected by emission or anomalous absorption, it 
is crucial to base spectra l type measurements o n several 
indice s (see discussio n in Herna ndez et alJ 2004) . Herbid 
(fl977l). IHerbig et all (fl986l) , iRvdgren et al.l (| 19761) , and 
ICohen fe Kuhil (|1979l) analyze spectral ranges of 5850- 
6700 A, 5800-6800 A, 3500-7000 A, and 4270-6710 A, 
respectively, with a resolution similar to that of our ob- 
servations, but with fewer spectral indices. The medium 
resolution (~2 A) spectra by lHartigan et a l. (1994) cover 
a much smaller range in wavelength (5700-6900 A and 
5860-6890 A) that includes few spectral indices. We in- 
clude more indices and these are based on a library of 
standards obtained with the same instrument and con- 
figuration. 

In Figures [TJ [2j and [3] we illustrate the accuracy of the 
SPTCLASS fitting procedure by comparing our targets 
to spectral type standards. In the late-type stars, the 
molecular bands of TiO and CaH were used in classify- 
ing the targets (Figures Q] and [2]) . In G-type stars, the 
metallic lines play an important role in spectral classifi- 
cation and the features used to classify UX Tau A include 
the Ca I, Mg I, Fe I, and Ca I lines (Figure EJ. The TiO 
and CaH bands in LkCa 14 (Figure^]) are in better agree- 
ment with a K5/K6 star. In GM Aur and LkCa 15, these 
bands are similar to those in a K5/K6 star and K4/K3 
star, respectively, rather than earlier (Figure [5]) or later 
(Figure Q]) type stars. The lines in UX Tau A, particu- 
larly Ca I and Na I, are weaker than those in a later type 
star (e.g. Kl) and resemble the depths of the lines in the 
G8 star but not the earlier type stars (Figure [3j) . Since 
there are no G6 or G7 objects in the library of standards 
used for our classification we cannot exclude it is a G6 or 



G7 star and so the uncertainty in UX Tau A's G8 spectral 
type classification is ±2.0. We note that our SpeX spec- 
tra are consistent with the spectral types derived with 
FAST. 

4.2. Near-Infrared Excess Emission 

To estimate the amount of veiling present in our sam- 
ple, we compared the line strengths in the target to those 
of a standard star of the same spectral type (5 I4.2J) . The 
spectra of corresponding template stars were ob tained 
from the IRTF Sp ectral Library (jCushing et alJ 120051 : 
iRavner et alJ [2009 j___. We first normalized the target 
and standard stars' emission at 2.2 /im and then arti- 
ficially veiled the standard spectrum by adding a flat 
continuum of F excess — r K x F standard where r K is the 
veiling factor (—F excess /F^). The standard spectrum was 
veiled until the depths of the absorption lines matched 
those of the target spectrum between 2.20 - 2.28 /xm. 
We normalized the standard spectrum in the K-band to 
the colors of a correspondin g stellar photosphere from 
iKenvon fe Hartmannl (|1995l ) scaled to the observed J- 
band flux. We then scaled the target by 1 + tk and 
subtracted the standard spectrum to extract the shape 
of near-infrared excess emission. When excess emission 
is present, we fit the continuum with a disk model and 
derive the best-fit parameters. 

4.2.1. LkCa U 

To test our method of measuring the veiling and 
excess emission of T Tauri stars (TTS ), we analyzed 
the spectrum of LkCa 14, a single (|White fe Ghezl 
120011), non-accreting diskle ss T Tauri star in Taurus 
( Kenvon fe Hartmannl Il995). The expectation is that 
this object should have no K-band veiling since it 
does not exhibit excess e mission at near- and mid- IR 
wavelength s in i ts SED (jKenyon fe Hartmannl |l995). 
iKarr et ail (|2010l ) report a tentative detection of dust 
around LkCa 14 based on H-band scattered light im- 
ages, ho wever, LkCa 14's emi ssion is photospheric out to 
~24 /urn (jLuhman et al.ll2010f ). making it suitable for the 
purposes of this work. 

We measure the amount of veiling in LkCa 14 by com- 
paring its K-band spectrum to that of a K5 V spectral 
standard from the IRTF spectral library (HD36003). We 
find that LkCa 14 has no veiling (r K =0±0.2) from 2.20 - 
2.28 fim. In order to extract any excess emission between 
2-5 ^m, if present, we first deredden the spectra. The 
standards from the IRTF library have already been cor- 
rected for reddening. For LkCa 14, we measure a visua l 
ext inction (Ay) of b y matching V jHjg et al.l [2000), 
R (jNortonet all 120071) . and 2MA SS photometry to the 
color s of a K5 stellar photosphere (jKenvon fe Hartmann 
1995) which is scaled to the J-band photometry. We sub- 
tract the standard from the target and find that there 
is no significant excess emission at longer wavelengths 
(Figure S]). There is a small mismatch between the two 
spectra from ^3.5-4.2 /xm which is most likely due to un- 
certainties in the measured spectral slopes of the target. 

We conclude that, as expected, LkCa 14 has no veiling 
and no excess emission. This indicates that for the pur- 
poses of this work dwarfs in the IRTF online library are 
suitable templates. 

13 http://irtfweb.ifa.hawaii.edu/~spex/IRTF_Spectral_Library/ 
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4.2.2. LkCa 15 



lEspaillat et al.l (|2008al ) measured a veiling factor at 
-2.2 /xm of 0.3±0.2 for LkCa 15. Using this veil- 
ing factor, they extracted the excess emission spectrum 
from the SpeX data and fit it with a simple black- 
body at a temperature of 1600 K, corresponding to the 
behavior expected from the wall of an inner optically 
thick disk. This blackbody fit, when taken in conjunc- 
tion with the dust clearing infe rred from the broad- 
band SED (Espaillat ct al. 2007a) and millimeter imag- 
ing (jPietu et al.ll2006l ). made LkCa 15 the first confirmed 
gapped disk around a TTS. 

Here we present new SpeX data from —0.8-5 /xm for 
LkCa 15 and analyze the near-IR excess emission of this 
ob ject in more deta i l by fi tting it with the disk wall model 
of lD'Alessioet al.l ([20051 ). We compare LkCa 15 a K3 
V spectral standard (HD219134) from the IRTF library 
and meas ure an ik of 0.6±0.2. While this is in agree- 
ment with Espailla t et al.l (|2008a|) within the uncertain- 
ties, it may also reflect infrar ed variability as has be en ob- 
served in transitional disks (Muzcrol le et al.l 120090 . We 
scaled the spectrum of LkCa 15 relative to the tem- 
plate according to this veiling factor and dere ddened the 
LkCa 15 spectrum using the iMathisI (J1990D deredden- 
ing law until it agreed with the template spectrum at 
1 fim. In this manner, we measure an extinction of 1.3. 
Extinction measurements based on matching V- , R- , 
I-band, and 2MAS S photometry to photospheric colors 
from Kcnyon & Hartmann (1995) derive Ay=1.7. If we 
scale LkCa 15's spectrum using r#-=0.4, which is within 
the veiling measurement uncertainty, we measure a visual 
extinction of 1.3 and so the difference in Ay is within the 
uncertainties of the veiling determination. Alternatively, 
this may be a reflection of the fact that in some cases one 
can derive higher visual extinc tions using color indice s at 
longer wavelengths than V-R (|Guilbring et al.l [1 9 981) . 

After subtracting the template spectrum from the 
target spectrum to extract the excess above the stel- 
lar photosphere, we fit the residual excess continuum 
with a disk wall model (Figure [5]). For the inner wall 
model of LkCa 15, we adopt the stellar prop erties listed 
in Ta ble [3] which were obtained with the iSiess et all 
(2000) evolutionary tracks. The mass accretion rate of 
Lk Ca 15 was estimated f rom its U-band excess follow- 
ing [GuIlbrmiliEjjI] (J1998I ). We adopt a minimum grain 
size of 0.005 /x m as is typically as sumed for the inter- 
stellar medium fMath i"s"et al.lll977l) and vary the max- 
imum grain size and the temperature of the optically 
thin wall atmosphere {T l wall ) to achieve the best-fit to 
the SpeX excess emission. We test wall temperatures be- 
tween 1000-2000 K, corresponding to the range of dust 
sublimation temperatures found for a large sample of 
classical T Tauri stars (CTTS) and Herbig Ae/Be stars 
(|Monnier fc Millan-Gabetl [2002), and maximum grain 
sizes of 0.25 /im, 1 /im, and 5 /xm. LkCa 15 has significant 
10 /xm silicate emission, indicative of small grains, and 
therefore we do not test larger grain sizes. The best-fit 
wall to the excess emission has a temperature of 1400 K, 
a height of 0.017 AU, and a maximum grain size of 1 /xm 
which, according to Eq. [51 corresponds to a wall located 
—0.12 AU from the central star (Figure [5]). This detailed 
model fit is consistent with the sim ple 1600 K blackbody 
fit found in earlier work (jEspaillat et al.l l2008al ) . The 



difference in temperature follows from the additional pa- 
rameters included in the wall model, namely the effect of 
the dust properties and the inclusion of an optically thin 
wall atmosphere. Our detailed fit to the excess contin- 
uum is con sistent a veiling measurement at —4.8 /xm by 
ISalvk et al.l ([2009T) (see Figure 15) . 

4.2.3. UX Tau A 

UX Tau A was first mo deled as a pre-transitional disk 
by lEspaillat et al. (2007a). To discern the nature of the 
innermost disk of this object, here we compare UX Tau A 
to a standard star of G8 V (HD101501) from the IRTF 
library and find that the absorption lines in the K-band 
spectrum of UX Tau A are weaker than those in the 
spectrum of the standard dwarf between 2.20-2.28 /xm. 
We measure an yk of 0.4±0.2 which is consistent with 
2MASS K-band photometry, from which we infer an tk 
of 0.3. 

We scale UX Tau A's SpeX spectrum relative to 
the template according to the derived veiling factor 
and deredden the target spectrum until it agrees with 
the template spectrum at 1 /xm. We measure an ex- 
tinction of 2.3 using this method, however, we de- 
rive an Ay of 1.8 for UX Tau A by comparing V- 
, R-, I-band (jKenvon fc Hartmann! 1 19951 ) and 2MASS 
hotometry to a standard G8 photosphere's colors 
Kenvon fc Hartmann! IT995). This difference in Ay is 
within the uncertainties of the veiling determination. 

After extracting the near-infrared excess of UX Tau A 
(Figure O, we fit it with the wall model described in 
§0 We use the stellar parameters (M*, R*, L») listed 
in Table [5] which were derived from the HR diagram and 
the ISiess et al.l (2000) evolutionary tracks, and we adopt 
an inclination for UX Tau A of 60°. As with LkCa 15, 
we adopt a minimum grain size of 0.005 /xm and vary 
the temperature of the wall between 1000-2000 K. Since 
there is no 10 /xm silicate feature in UX Tau A, indicating 
a lack of small grains, we adopt a maximum grain size of 
10 /xm. We cannot discriminate between wall fits using 
a maximum grain size of 10 /xm, 100 /xm, or 1000 /xm, 
however, since the wall's parameters (z^ aZi , R° all ) do not 
vary significantly between the three grain sizes for the 
purposes of this work, we adopt the smallest maximum 
grain size that can reproduce the SpeX excess spectrum. 
We find that the best-fit wall has a surface temperature 
of —1550 K, is located 0.15 AU from the star, and has a 
height of 0.009 AU (Figure |S}. Our results are consistent 
with other veil i ng me asurements found in the literature. 
Edwards et al.l (120061) measure a veiling of at 1 /xm and 
Salvk et al . (2009)'s veiling measurement at —4.8 /xm is 
consistent with our analysis within the veiling measure- 
ment uncertainties (see Figure [5]). 

4.2.4. Rox 44 

Rox 44, located in the Ophiuchus cloud, was identified 
as a pre-transitional disk candidate based on its strong 10 
/xm silicate emission f eature which could not be explained 
by fu ll disk models (jFurlan et al.l 120091: IMcClure et"aTl 
2010). Rox 44's Spitzer IRS spectrum is also strikingly 
similar in shape to that of LkCa 15. To investigate the 
nature of the inner disk of Rox 44, we obtained a near- 
IR SpeX spectrum. When compared to a K3 V standard 
(HD219134) from the IRTF library, the absorption lines 
in the K-band spectrum of Rox 44 are weaker between 
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2.20 - 2.28 /im and we measure an xk of 0.8±0.2, con- 
sistent with an tk of 0.7 inferred from 2MASS K-band 
photometry. 

We scale the spectrum of Rox 44 relative to the tem- 
plate according to this veiling factor and deredden its 
spectrum until it agrees with that of the template spec- 
trum at 1 /im (Figure [5]). We measure a visual extinction 
of 2.2, in a greement with previous extinction measure- 
ments from lBouvier fc A ppcnzcllcr (1992)0 The excess 
emission above the stellar photosphere can be fit with 
a wall model (Figure [5]) using the stellar parameters in 
Tableland adopting a minimum grain size of 0.005 /im. 
We tested maximum grain sizes of 0.25 /im, 1 /an, and 
5 /im and temperatures between 1000-2000 K. The best- 
fit wall to the near-infrared excess continuum of Rox 44 
has a height of 0.034 AU, a maximum grain size of 1 /im, 
a temperature of 1200 K, and is located 0.25 AU from 
the central star (Figure [5]) . 

4.2.5. GM Aur 

We now turn to the transitional objects in our sam- 
ple. The transitional disk around the sin gle star GM Aur 
(jPott et all [2010t IWhite fc Ghe3 l200l A. Kraus, pri- 
vate communication ) has an inner hole of ~ 20 AU 
([Hughes et al.l l2009f ) wh ich contains some sub-micron- 
sized optically thin dust (jCalvet et al.ll2005h . In our veil- 
ing analysis, we compare GM Aur to a K5V standard 
(HD36003) and find that GM Aur's K-band absorption 
lines are slightly veiled (rK=0.1±0.2) and the object ex- 
hibits a small excess above the stellar photosphere (Fig- 
ure [5j) . This veiling measurement is consistent with a 
veiling factor of 0.2 inferred from 2MASS K-band pho- 
tometry. 

After scaling the GM Aur spectrum according to 
its veiling factor, we derive an extincti on (Ay) of 1.2 
by dereddening its spectrum with the Mathis (1990) 
law until it matches the template spectrum at 1 /im. 
However, as in the case of LkCa 15 and UX Tau A, 
this is higher than the extinctio n of 0.8 measured by 
matc hing V, R, and I photometry (jKenvon fc Hartmannl 
19951). and 2MASS fluxes to photospheric colors from 



Kenvon fc Hartmannl (J1995I ). 



The small excess seen in GM Aur's SpeX spectrum is 
in line with the excess emission ab ove the photospher e 
seen in its SEP fr om ~ 3-8 /im (jCalvet et al.l I2005T ). 
iCalvet et al.l (J2005D showed that this emission comes 
from optically thin sub-micron-sized dust located within 
the inner 5 AU of GM Aur's 20 AU hole. We can fit the 
SpeX excess of GM Aur with an optically thin dust model 
Figur e [5]) supporting the conclusions of ICalvet et al.1 
2005). The optically thin dust region lies within ~1 AU 
of the star and is composed of ^2xl0 -12 M© of dust 
made up of silicates (47%), organics (42%), and troilite 
(11%). The total emission of this optically thin region is 
scaled to the vertical optical depth at 10 /im (to ~ 0.021). 
We note that since the structure of the dust within the 
holes of transitional disks is not known at this time, our 
results should be taken as an approximation of how much 



14 This ex t inctio n differs from the Ay=3.3 reported by 
McClurc et al. (2010) based on near-IR colors. In this work we 
adopt an Ay that agrees with the optical photometry. More opti- 
cal photometry is needed to better determine the extinction of this 
object. 



dust is needed within the hole to reproduce the SED, not 
as a detailed fit from which we can derive the spatial dis- 
tribution of the optically thin dust near the star. High- 
resolution near-IR interferometry is needed to trace this 
component in more detail. 

Our results are in agreement with other veiling mea- 
surements at 1 nm (Edwa rds et al.l I2006D and 2.2 /im 
(jFolha fc Emersor] Fl999'). We note that the excess emis- 
sion of GM Aur derived here does not agree w i th th e 
J-band veiling measurement of iFolha fc Emerson! (J1999). 
However, these authors use a later spectral type standard 
(K7) in their analysis which could have resulted in the 
higher veiling measurement since later type stars have 
deeper absorption lines. Alternatively, this may be a re- 
sult of variability as i s seen in other transitional disks 
(jMuzerolle et al.li2009l ). 

4.2.6. DM Tau 



DM Tau is a single star (jPott et al.l 120101: 

IWhite fc Ghezl 120011 A. Kraus, private communica- 
tion) surrounded by a transitional disk wi th a ~3 AU 
hole that is empty of small dust grains (jCalvet et al.l 
2005). The absorption lines of DM Tau are not veiled 
relative to a Ml V standard (HD42581; Figure EJ) and 
we measure an r#- of 0±0.2 between 2.20 - 2.28 /im in 
the SpeX spectrum. 

We derive an extinction of (Ay) of 0.7 for DM Tau by 
fittin g V-, R-, I-band photometry (jKenvon fc Hartmannl 
I1995D and 2MASS colors to an M l stellar photosphere 
from iKenvon fc Hartmannl (J1995I ). Th e spect r um is 
dereddened with this extinction and the Mathis (1990) 
dereddening law. We note that we did not obtain a low- 
resolution SpeX spectrum at shorter wavelengths for this 
target. Subtraction of the standard from the target re- 
veals that DM Tau exhibits no significant near-IR excess 
emission at 2-5 /im (Figure [6]), in agreement with its lack 
of emission above the photosphere in the K-band (see 
SED in Calvet et al. 2005 or Furlan et al. 2006) as well 
as a veiling measurement at ~4.8 /im that i s consistent 
with no excess emission at this wavelength (jSalvk et al.1 
12009( 1. As in the case of LkCa 14, DM Tau shows a 
small mismatch when compared to the standard, which 
may be an indication of spectral slope uncertainties. The 
lack of veiling and near-IR excess emission in DM Tau 
is consistent with the fact that this object has no excess 
emission above the photosp here from the K-band up to 
- 8 /im (jCalvet et al.ll2005f i. indicating that DM Tau has 
relatively no small dust within its hole. 

4.3. Broad-Band SED Modeling of Pre- Transitional 

Disks 

Here we model the broad-band SEDs of the pre- 
transitional disks of LkCa 15, UX Tau A, and Rox 44 
with the best-fit inner wall model from § 14. 2| includ- 
ing the shadowing of this inner wall on the outer disk. 
LkCa 15 and UX Tau A have been modeled previously 
(jEspaillat et al.|[2007al) . but here we incorporate their 
ne wly d erived stella r properties (S I4.1J) . and inner walls 
(§ 14X21 and g l4~2~3|) . and a finite star. 

4.3.1. LkCa 15 

For LkCa 15, we measure a 58 AU gap that contains 
some small optically thin dust. This is compatible with 
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millimeter interferometric imaging bv lPietu et al.l (|2006T ) 
which finds a cavity of 46 AU in this disk, especially 
when one takes into account the uncertainties introduced 
by dust opacity one assumes (see § I5.2J) . Figure [7] con- 
tains the broad-band SED of LkCa 15. Observations 
include ground-based optic al photometry (open circles; 
iKenvon fc Hartmann|[l99 5-). and J, H, K (filled circles; 
2MASS), SpeX (red s olid line; this work) Spitzer IRAC 
(blue solid trian gles; ILuhman et al.l [2010fl. Spitzer IRS 
(blu e solid line: IFurlan et al.l l2006f ) . IRAS (open trian- 
gles; l^feaverjklJoncs 1992), and millimeter (filled pen- 
tagon from Andrews & Williams 2005; open pentagons 
from Espaillat et al. 2007) data. The SpeX spectrum 
used in this paper is not flux calib rated; here we scale it 
to th e adopted stellar photosphere fKen von fc Hartmannl 
I1995T ) using the veiling derived in § 14.2.21 

In Figure [3 we use the best-fit wall model found in 
§ 14.2.21 We note that while this wall is a good fit to 
the SpeX spectrum, it over-predicts the Spitzer IRAC 
and IRS data. This can be explained by the uncertain- 
ties inherent in our veiling measurement or the fact that 
infrared variability in these types of disks is common (Es- 
paillat et al., in preparation: iMuzerolle et al. 2009). The 
10 /jm silicate emission can be fit by an optically thin 
dust model which contains — 2xl0~ n M© of optically 
thin dust within the inner 4 AU of the gap. This re- 
gion is composed of 79% silicates, 12% organics, and 9% 
troilite and the total emission is scaled to the vertical 
optical depth at 10 /zm (to ~ 0.009). The outer disk 
has a=0.0005 and e=0.001 and is inwardly truncated at 
about 58 AU. The inner wall dominates the near-IR emis- 
sion and we constrain the inner disk to <.19 AU in radius, 
assuming that the inner disk has the same a and e as the 
outer disk. Model parameters are listed i n Table El 

Our results for LkCa 15 differ from Espaillat et al. 
(l2007aT). In fir s t mod eling the outer disk of LkCa 15, 
lEspaillat et al.l (|2007al) adopted a spectra l type of K5 
(T*=4350 K) for this object following iHerbig et all 
(1986). However, the spectral type we derived in § 14.11 
for LkCa 15 (K3) corresponds to a higher temperature 
(T*=4730 K) than previously adopted. The largest dif- 
fe rence between the disk model presented here and that 
in Espaill at et al.l (|2007a|) is the location of the outer disk 
wall. Assuming a spectral type of K5 the disk wall was 
located at 46 AU while for a K3 star the wall is located 
at 58 AU. This difference is expected since the location 
of the wall, given by Eq. [21 depends on L* and L acc and 
these parameters are now higher with respect to the pre- 
vious model. 

Based on the equations presented in the Appendix, the 
inner wall of LkCa 15 would cast an umbra of about 7 AU 
on the outer wall (measured from the midplane) in the 
case that the star is taken to be a point source. If the star 
is taken to be a finite source, LkCa 15's outer wall has 
an umbra of —4 AU and the upper edge of the penumbra 
is located at -9 AU (Eq.EH&EU)- The result for the 
solid angle (Eq. IA6[) of the star seen by an observer at 
different heights of the outer wall of LkCa 15 is shown in 
Figure [8] Taking into account the shadowing of a finite 
star, the outer wall of LkCa 15 is 12.9 AU high. This is 
2.6 H, where H is the gas scale height, and is consistent 
within 15% of the independently calculated height of the 
outer disk (11.7 AU) at -58 AU. Therefore, while the 
outer wall is partially shadowed by the inner wall, we 



conclude that the outer disk behind the wall is not com- 
pletely shadowed by the inner disk. The height of outer 
wall in this wo r k is hi gher than previously reported by 
lEspaillat et al.1 ()2007al ). The z wa u in that paper repre- 
sents the part of the wall that is illuminated when the 
star is taken to be a point source, but was mistakenly 
defined as the height of the wall above the midplane. 

4.3.2. UX Tau A 

For UX Tau A, we find a 71 AU gap that is relatively 
devoid of small dust grains, as indicated by its lack of a 
10 /im silicate feature (Figure [9]). Observations shown in 
Figure E] i nclude ground- based optic al photometry (open 
circles: IKenvon fc Hartmannl 119951) . J, H, K photome- 
try (filled circles; 2MASS), a SpeX spectrum (solid red 
line; this work), Spitzer IRA C photometry (blue solid 
triangles; ILuhman et al.ll2010[ ) , a Sp itzer IRS spectrum 
(blue solid line; IFurlan et al.l 120 06) , IRAS photometry 
(open triangles; iWeaver fc Joneslll992j). and millimeter 
data ( pentagons; lAndrews fc Wil liams 2005|). Our SpeX 
spectrum is not flux-calibrated and we scale it to the 
adopted stellar photosphere using the veiling derived in 

§0231 

We use the best-fit wall model found in I 



1 l4T3l in fitting 
the SED emission presented in Figure \§\ We find that 
the outer disk has a=0.004 and e=0.001 and is inwardly 
truncated at about 71 AU. The inner wall dominates 
the near-IR emission and we constrain the inner disk, 
assuming it has the same settling and viscosity as the 
outer disk, to <.21 AU in radius. Model parameters are 
listed in Tabled 

The outer wall of UX Tau A would have an umbra of 
about 4 AU in the case that the star is taken to be a 
point source. If the star is taken to be a finite source, 
the size of the umbra is —0.4 AU and the upper limit 
of the penumbra is —8 AU. The result for the solid an- 
gle seen by an observer at different heights of the outer 
wall of UX Tau A is shown in Figure [5] The outer wall 
of UX Tau A is 13.8 AU (2.0 H) high and is consistent 
within 15% of the height of the outer disk (12.5 AU) at 
~71 AU, which is calculated independently. The fully 
illuminated part of the outer wall is at a temperature 
where ice has sublimated (110 K), however, in the par- 
tially shadowed portions of the wall the temperature is 
lower. In order to make the composition of the wall con- 
sistent, we remove ice from the partially shadowed por- 
tions of t he wall. We note that UX Tau A has crystalline 
silicates ([Espaillat et al.ll2C)07aI ) , which we do not include 
here and would contribute to the IRS spectrum between 
~20 - 35 /im. 

As in the case of LkCa 15, our results for UX Tau A dif- 
fer from lEspaillat et al.l (|2007af ) because we are adopting 
a different spectral type in this paper. Es paillat et al.1 
(|2007a| )'s model was based on a s pect ral type of K2 
T*=4 900 K) followingEerbii (fl977l) and lCohen fc"Kuhl 
19791 ). The spectral type we derived in § 14.11 for 
UX Tau A (G8) corresponds to a significantly higher 
temperature (T*=5520 K) than previously adopted. The 
difference in R° all is expected following Eq. [2] 

4.3.3. Rox 44 

Rox 44 has a 36 AU gap that contains some ISM- 
sized optically thin dust. The size of this gap is con- 
sistent with the 33 AU cavity observed with the SMA 
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(|Andrews et al.ll2009D . Observations shown in Figure _ 
inclu de ground-based optical (open circles; iHerbst et al 
Il99l . J, H, K (filled circles; 2MASS), Spe X (solid red 
line; this work), Spi tzer IRAC and MIPS (jEvans et al.1 
l2009HPadgett et all 120081 blue solid tria ngles), Spitzer 
IRS (blue solid line; [Fu rlan ct al. 2006) , and millime- 
ter (p entagons; lAndrews et al.l 120091 : iNuernberger et al.l 
1998) data. The SpeX spectrum is scaled to the adopte d 
stellar photosphere using the veiling derived in §14.2.41 

We use the best-fit wall model found in § 14.2.41 in 
the SED model fit presented in Figure [TO] Rox 44 has 
~2xl0 -11 M of optically thin dust within 2 AU that 
contributes to the 10 fim silicate feature and that is com- 
posed of 79% silicates, 12% organics, and 9% troilite. 
The total emission of this region is scaled to the vertical 
optical depth at 10 /im, To ~ 0.036. The outer disk has 
a=0.006 and e=0.01 and is inwardly truncated at about 
36 AU. The inner wall dominates the near-IR emission 
and we constrain the inner disk to <.4 AU in radius. 
Model parameters are listed in Table [3) 

In the case that the star is taken to be a point source, 
the outer wall of Rox 44 would have an umbra of ~5 AU. 
For a finite-source star, the size of the umbra is ~4 AU 
and the upper edge of the penumbra is located at ~6 AU 
(see Figure |SJ|. Taking into consideration this shad- 
owing of the outer disk, the outer wall of Rox 44 has 
z w;~9.9 AU (3.5 H) and this height is consistent within 
~35% with the height of the outer disk at ~36 AU 
(7.3 AU). We note that far- infrared variability is likely 
in this object given that the Spitzer IRS and MIPS data 
were taken at different epochs and are not in agreement, 
even when taking the uncertainties of the observations 
into consideration. If the far-infrared flux of Rox 44 was 
less than is seen in Figure [10] at the time when the SpeX 
observations were taken, then we would be overestimat- 
ing the height of the outer wall. Simultaneous observa- 
tions of this object in the near- and far-infrared would 
be useful in testing this. Additional far-infrared and 
sub-millimeter data would also help further constrain the 
properties of the outer disk. 

5. DISCUSSION 
5.1. Pre- Transitional Disk Structure 

It is now accepted that full disks have a sharp transi- 
tion at the dust destruction radius, inside of which the 
temperature is too high for dust to exist. This tran- 
sition appears as a "wall" that is frontally illuminated 
by the star, since for typical m ass accretion rates th e 
inner disk gas i s optic ally thin ([Muzerolle et al J 120041 ). 
iMuzerolle et al.l (|2003| ) found that nine T Tauri stars 
surrounded by full disks in Taurus had SpeX near-IR 
excess emission which could be fit by single temperature 
blackbodies with temperatures that fell within the range 
of dust sublimation te mperatures (1000-2000 K ) foun d 
for a larger sample by iMonnier fc Millan-Gabetl (|2002t ). 
IMuzerolle et al.l (|2003l )'s fits are evidence that there is 
optically thick material located at the dust destruction 
radius in these full disks, most likely from the inner wall 
of the disk which emits primarily in the near-IR. 

Pre-transitional disks exhibit the same behavior as full 
disks in the near-IR. In this work, we can fit the continua 
between 1-5 /zm of the pre-transitional disks in our sam- 
ple with an inner disk wall which is located at the dust 



destruction temperature. This conclusion is consistent 
with independent veiling measurements found in the lit- 
erature (see § I4.2[) as well as near-IR interfero metric mea- 
surem ents. Using the Keck interferometer, iPott et aLl 
(2010) find that the K-band emission of LkCa 15 and 
UX Tau A originates from very small radii. This is in 
good agreement with what we report in this wo r k (Ta - 
ble [3]) if one takes into account that IPott et aT] ()2010f ) 
used simple models assuming a face-on disk and the de- 
rived radii depend on the dust opacities that are adopted. 
Independent SED modeling also s upports that the inne r 
disk of LkCa 15 is optically t hick JMulders et aTll2010D . 
We note that iMulders et~all ()2010D claim that the opti- 
cally thick inner disk of LkCa 15 extends out to 1 AU 
based on fitting broad-band photometry and assuming 
that the shadowing in the disk is due to a p oint source. 
This r adius is inconsistent with the results of IPott et al.1 
(2010) and we find that an inner disk of this size does 
not fit the Spitzer IRS spectrum. 

While pre-transitional disks have the same behavior as 
full disks in the near-IR, at longer wavelengths in their 
SEDs and in millimeter images, they display evidence 
for large dust clearings within the disk. In the Spitzer 
IRS spectra of LkCa 15, UX Tau A, and Rox 44, we see 
a pronounced decrease in the emission at ^20/im which 
indicates a drop in the opacity of the disk, most likely 
due to the removal of dust. This is supported by millime- 
ter in terferometric observations of LkCa 15 (jPietu et al.1 
I2006T ). UX T au A (S. Andrews, pr ivate communication), 
and Rox 44 ((Andrews et al.1 120091 ) which see cavities in 
the images and sharp drops in the visibilities, evidence 
that dust in the disk has been removed. In conjunction 
with the evidence for an optically thick inner disk, this 
leads one to the conclusion that pre-transitional disks 
have gaps in their disks. 

These gaps in pre-transitional disks are clearly differ- 
ent from the inner holes detected in transitional disks (see 
Figure ITT1 for a schematic). DM Tau has an inner hole of 
3 AU which is relatively clear of small dust as indicated 
by the lac k of near-infrared emission in its Spitzer IRS 
spectrum (jCalvet et al. 2005(). This is supported by the 
lack of veiling in the K-band as well as photospheric near- 
infrared emission, resembling what is seen in the diskless 
LkCa 14, indicating that the inner disk of DM Tau is de- 
void of small dust. GM Aur is another transitional disk 
studied in this pa per and its SED indicates that it has 
a ~20 AU hole (jCalvet et al.ll2005 | ) which ha s been con- 
firmed in the millimeter (|Hughes et "all 120091 ). GM Aur 
has a small excess above the photosphere in the near- 
infrared, as measured by its Spitzer IRS spectrum, which 
has been interpreted as originating from optically thin 
dust located within the disk hole (Cal vet et al.l 12005). 
In this work we demonstrate that the small, flat near- 
IR excess continuum extracted from the SpeX spectrum 
can be fit by a model with ~10~ 12 M Q of optically thin 
dust within 1 AU of the star. About 85% of the emis- 
sion in the K-band originates in the inner 0.5 AU of the 
disk, consistent with ncar-intcrfcromctric measurements 
(|Pott et al.ll2010l ). 

An alternative model for t he inner disk of pre- 
transitional disks propose d by Espailla t et all (|2007al) 
and further investigated by Mu lders et al.l (|2010l) is that 
the inner disks of these objects are optically thin. How- 
ever, GM Aur's near-IR emission is strikingly differ- 
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ent from that seen in our pre-transitional disk sample. 
GM Aur's near-IR excess is relatively flat, reflecting a 
range of different temperatures. On the other hand, the 
near-IR excess emission of pre-transitional disks resem- 
bles a single temperature blackbody. Optically thin dust 
spread out over several radii will not appear as a single 
temperature blackbody and if one confines this optically 
thin dust to very small radii, one would need a shell of 
dust engulfing the star ([Mulders et al.l [2010) for which 
there is yet no supporting evidence. 

Since the inner disk in pre-transitional disks is opti- 
cally thick, it will cast a shadow on the outer disk. The 
shadow cast by the inner wall on the outer wall has sig- 
nificant consequences on the sizes of gaps that we can 
detect with Spitzer. The majority of the emission traced 
by IRS comes from within the inner 10 AU of the disk 
and ov er 80% of the emission at 10 /xm comes from within 
I AU (D'Ales sio et al.ll2006t lEspailladl2009T ). Therefore, 
the Spitzer IRS instrument will be most sensitive to clear- 
ings in which some of the dust located at radii <1 AU has 
been removed. To test the limits of the smallest gap IRS 
could detect around a K-type star, we simulated gaps 
in a disk around a 0.5 Mq star with a mass accretion 
rate of 10 -8 M Q yr _1 , an e of 0.01, and an inclination 
of 60°. In this disk, we find that the smallest gap that 
will cause a noticeable "dip" in the Spitzer spectrum will 
be ~4 AU. As is seen in Figure [T2l the SED of such 
a gapped disk begins to differ from that of a full disk 
beyond 10 /xm. This is because, relative to a full disk, 
there is extra emission from the wall of the outer disk. 
Since the height of the disk increases with radius, the 
inner disk can only extend from the dust destruction ra- 
dius (0.12 AU) out to about 0.3 AU before it significantly 
shadows the outer wall, thereby dimming the emission of 
the outer wall and making it indistinguishable from a full 
disk. Furthermore, IRS cannot easily detect gaps whose 
inner boundary is o utside of 1 AU (e.g. a gap spanning 
5-10 AU in the disk; lEspaillati[2009h and it will be up to 
ALMA to detect such gaps. 

There are some caveats to keep in mind regarding the 
above estimate of the IRS gap detection limit. First, we 
are assuming that the inner disk has the same amount of 
dust settling as the outer disk. However, if the inner disk 
is more settled it will have a lower surface and will not ob- 
scure the outer wall as much, and so the inner disk could 
extend out to further radii. Also, our simulations do not 
include optically thin dust within the gap and we note 
that we did not include any possible shadowing of dust 
in the optically thin regions inside the gaps of LkCa 15 
and Rox 44. Infrared interferometry is crucial to explore 
how the spatial distribution of this component is shaped 
as it filters into the gap from the outer disk (iRice et al.1 
120061) or how it is i nfluenced by possible planets in the 
gap (JLubow fe D'Angeldl200l . 



5.2. Dust Opacities 

The opacity of the disk is controlled by dust and in any 
sophisticated disk model the largest uncertainty lies in 
the adopted dust opacities. The dust opacity sets the lo- 
cation of the disk surface, i.e. where the optical depth to 
the stellar radiation reaches ~ 1 , and in turn the amount 
of disk flaring. The shape of the disk surface dictates the 
fraction of the stellar radiation that is captured by the 
disk, and as a result, the heating of the disk and the 



ensuing disk emission (jCalvet et al.lll99ll) . 

The SED models derived in this work are obviously 
tied to the dust opacity that is assumed in § [3J In par- 
ticular, R° aH depends on the adopted opacity follow- 
ing Eq. [2l We find t hat b y adopting the dust opac- 
ities of iDraine fc Led (|1984[ ). the derived R° Q ;/ of our 
targets increases by ~10 AU. This occurs because the 
opacity to the stellar radiation (k s ), which peaks around 
1 ^m, is higher for the dust composition of lDraine fc Led 
(1984) than that adop ted in this work (Figure fT3|) . 
For I Pollack et al.l (|1994t )'s opacities, R° aW decreases by 
~10 AU. In this case, the opacity to the disk's local emis- 
sion (Kd), which peaks ^25 /im, is much higher than ours. 

The major diff erences between the opaci ties we adopt 
here a nd those of lDraine fe Led (|1984f) and lPollack et all 
(1994) are the assumed carbon component and the ice 
abundance. While most models assume some carbon 
component in their dust composition, there is yet no 
general consensus o n the species of the carbon grains. 
IDraine fc Led (|1984fl a ssume that the carbo n component 
is due to graphite and iPollack et al.l (|1994f ) assume it is 
due to organics. Further study is needed to better deter- 
mine the carbon component in accretion disks. The Her- 
schel mission is an excellen t opportunity t o lear n more 
about water ice in disks. IPollack et al.l (|1994f) adopt 
a dust-to-gas mass ra tio for water ice of 0.0056 while 
iD'Alessio et al.l (|2006l ) find that the features produced 
by this value do not agree with observations of disks. 
Herschel can help settle this issue by detecting water 
lines in the far-infrared and allowing derivations of the 
abundance of water in disks. 

5.3. Implications of Pre- Transitional Disks on Dust 
Clearing Mechanisms 

Several mechanisms have been proposed to explain disk 
clearing. These include planets, grain growth, the mag- 
netorotational instability, and photoevaporation by the 
central star. The structure of pre-transitional disks - 
an inner disk, a gap, and an outer disk - is predicted 
by planet formation models. On the other hand, the 
three alternative clearing mechanisms struggle to make 
a gapped disk structure and, due to their predominantly 
inside-out clearing pattern, are instead more likely to 
form a central hole in a disk. 

Given the simple stipulation of a disk gap, planets 
emerge as the most likely dust clearing mechanism in 
pre-transitional disks. To further explore this possibility, 
one must ascertain if planet formation timescales are in 
line with the ages of these disks, if there is enough mass 
in these disks to have formed planets, and if planets can 
open the kind of gaps that we are detecting. Assuming 
that planets are indeed present within these disk gaps, 
one can attempt to place some constraints on theoretical 
planet formation models. 

It is feasible to form planets within 1 Myr, the age of 
our pre-transitional disk sample, via the two leading the- 
ories of planet formation: the gas instability model and 
the core nucleated accretion model. In the gas instabil- 
ity model, fragmentation due to gravitational ins tabili- 
ties in the dense dust midplane (|Goldreich fc W ard 1973) 
can form Jupiter-mass clumps in t h e outer disk within 
a few hundred ye ars (jBossI l2000t iMaver et all 120021 ; 
iDurisen et al.l 120071). According to the core nucleated 
accretion model ( Pollack et al.| [l996; Bodenhcimc r et al.l 
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2000) , dust grains will grow into planetesimals which will 
accrete other planetesimals and then form solid cores. 
These cores can become terrestrial planets or, if the 
cores are massive enough, gas giant cor es which can 
accrete the surrounding gas i n the disk ([Pollack et al.1 
1996; Bodcnhe imer et al.ll2000h . For several years it was 
thought that the formation time for giant planets via 
core accretion was longer than the lifetime of disks, how- 
ever, more recent simulations show that giant planets can 
form in the inner disk within ~1 Myr via two different 
mechanisms. One possible mechanis m which could accel - 
erate planet formation is migration. lAlibert et al.l (2005) 
found that a migrating core will not deplete its surround- 
ing feeding zone and will reach its cross-over mass, the 
point at which the planet enters runaway gas accretion, 
sooner than found in previous models. Alternatively, 
the time to enter runaway gas a ccretion can be red uced 
by adopting low dust opacities (Pollack et al.l 119961 ) . A 
Jupiter-mass planet wi th a core of 10 Earth -masses can 
form at 5 AU in 1 Myr ([Hubickvi et al.ll2005lh about half 
the time predicted by models using higher dust opacities. 

It is likely that the pre-transitional disks originally 
had enough mass within their presently cleared regions 
to have formed multiple giant planets in situ and that 
planet formation may still be underway. Compared 
to a full disk with the same model parameters (Ta- 
ble EJ, LkCa 15, UX Tau A, and Rox 44 are "missing" 
10 M j U piter, 40 M j U p iter , and 20 M Ju p lter , respectively, 
some fraction of which could have gone into making giant 
planets. There is still enough material in the outer disk 
of these objects (Table [3]) to potentially go into forming 
planets. 

The next issue to address concerns the nature of 
the signatures that forming planets will leave be- 
hind in young disks. Theories predict that newly 
forming planets should clear the material around 
themselves throu gh tidal disturbances leavi n g be- 
hind disk gaps ( Goldreich fc Tremainel 119801: iWardl 



19881: iRice et all 120031: iPaardekooper fc Mellemal 120041: 
Quill cn et al.l l2004t iVarniere et a l. 2006). To clear out 
the relatively large gaps we measure, it is likely that more 
than one planet is present and detailed simulations on the 
effects of multiple planets on the dust distribution are 
necessary. It is already known that when a stellar com- 
panion is loc ated in the inner disk, the dust is dynami- 
cally cleared ([Mathieu et al.l[l99ltlArtymowicz fc Lubowl 
1994). One example is the cavity in the inner disk of 
CoKu Tau/4 which i s most likely due to clearing by 
its stellar companion ([Ireland fc Krausll2008l) and whose 
SEP has been repr oduced by a circumbinary disk model 
(|Naeelet al.ll2010D . The object T54 in the Chamaeleon 
cloud is also a lik ely candidate fo r dust clearing by a 
stellar companion (Ki m et aLll2009l ). These observations 
are evidence that dynamical dust clearing by companions 
does occur; the salient question is whether the companion 
in pre-transitional disks is in the stellar-mass or planet- 
mass regime. Searches for companions in the disks of 
UX Tau A and LkCa 15 have revealed that they are single 
stars down to 0.35 AU (pott et aJJl2010b iWhite fc Ghezl 
1200 ll A. Kraus, private communication). To explain the 
location of the outer walls of UX Tau A and LkCa 15, 
a stel lar companion would have t o be located at larger 
radii ([Artvmowicz fc LubowllT9"94l ) , therefore we can con- 
clude that the truncation of the outer disk is not due to 



a stellar companion. Multiplicity studies of Rox 44 down 
to similar small radii do not yet exist. 

If planets are clearing the dust in pre-transitional disks, 
we can speculate on the possible masses of such planets 
and their most likely formation scenario. Simulations 
of disk clearing demonstrate that Neptune-mass planets 
(^0.05 Mj U pi ter ) can open gaps in the dust component of 
the disk while planets with 0.5 Mj UT>i t e r can cre ate a gap 
in the gas disk (Paardc kooper fc Mellemall2006|) . We can 
conjecture that the planets in the pre-transitional disks 
are between 0.05 Mj up i ter .-0.5 Mj up i ter based on the de- 
tection of a gap in their dust disks and evidence which 
suggests there are not substantial gaps in the gas disk. 
Since these disks are actively accreting, we can assume 
that gas travels from the outer disk to the inner disk 
and that therefore gas is located within the dust gap. 
In fact, in LkCa 15, UX Tau A, an d Rox 44, gas has 
been detected wi thin the dust gap ([Naiita et al.l 120031 : 
iSalvk et "aLl I2009H . These partially- filled gas gaps have 
consequences on Type II migration. According to classi- 
cal Type II migration, planets will migrate inward onto 
the star on the viscous timescale, followed by the outer 
disk edge, and th erefore we would not see gaps on these 
short timescales ([Lin fc Papaloizoul 119861 ). The gap in 
the gas disk in the above sce nario is deep and there is no 
gas flow across it. However, iCrida fc Morbidelhl ([2007D 
showed that inward planet migration can be slowed or 
stopped if the gap is not completely cleared of gas. In 
this case, the gas in the gap will exert a positive torque 
on the outer disk which diminishes the total torque that 
the planet feels from the outer disk. As a result the outer 
disk will not migrate the planet inward as quickly, open- 
ing the possibility that planets can sustain gaps in disks 
long enough for us to detect them. Based on the points 
above, it follows that the planets in pre-transitional disks 
were most likely formed via core accretion since gravita- 
tional instability tends to form planets of several Jupiter 
masses in the outermost disk, and moreover, these plan- 
ets would have needed to migrate to the inner disk to 
form the gaps we are detecting. 

The alternative disk clearing mechanisms proposed to 
date have difficulties explaining several observations and 
arc most likely not responsible for the type of clear- 
ing seen in pre-transitional disks. Grain growth simula- 
tions predict that dust will evolve on faster timescales 
in the inner parts of the disk, and eventually the 
small grains which contribute to the near-IR emission 
of the disk will grow causing a flux deficit in the SED 
(|Dullemond fc Domini]^ 120041 . 120051 ) . However, the gap 
edges and silicate features are too sharp to be explained 
by grain growth. In addition, grain growth is an inside- 
out clearing mechanism and does not account for a rem- 
nant inner disk. 

In the MRI clearing scenario, X-rays from the star can 
activate the MRI in the ionized inner wall of the disk 
which will lead material to accrete from the wall onto 
the star, creating a hole in the disk that grows outward 
([Chiang fc Murray-Clavi 120071 ). It predicts a stronger 
trend of gap size wit h X-ray luminosi ty and accretion 
rate than is observed ([Kim et al.1 120091 ) and requires an 
inner disk hole, not a gap, to already exist before it can 
take effect. 

Radiation from the central sta r can potentially photo- 
evaporate the surrounding disk (Hollcnb ach et al.lll994t 
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IClarke et alj|200lh . High energy photons from the stel- 
lar wind will impinge upon the upper disk layers. As 
the disk evolves viscously, the mass accretion rate de- 
creases with time, eventually reaching the mass loss rate 
in the photoevaporative wind. At this point, the photo- 
evaporative wind takes over and inward accretion onto 
the star will essentially stop. Once the isolated inner 
disk drains onto the star on the viscous timescale, the 
inner edge of the outer disk will be directly irradiated 
by the star and the hole w ill keep growing outward 
f Alexande r fc Armi tagc 2007J). The typical value of the 
mass loss rate in the photoevaporative wind was typ- 
ically taken to be ~4xl0~ 1 ° M w yr -1 based on EUV 
radiation alone ([Clarke et al.H200JT) . however, this value 
has recently jumped to ^10 -8 M^ yr -1 due to the inclu- 
sion of X-ray and UV photons ([Gorti fc Ho llcnbach 2009; 
lOwen et al.ll2010i ). However, if mass loss rates are as high 
as this value it poses significant problems. The average 
mass accretion rate me asured in TTS is ~10~ 8 M Q yr _1 
(jHartmann et al.|[l9 98). According to the evolutionary 
scenario outlined above, it follows that we should not 
see disks with mass accretion rat es lower than this value , 
in contrast to what is observed (Hartma nn et al.l ll998). 
Moreover, modeling of disk wind indicato rs (i.e. forbid- 
den line profiles) by IHartigan et all (J1995D find that the 
winds of TTS are typically ~0.01 of the mass accretion 
rate. In any case, photoevaporation opens short-lived 
gaps in the disk that would not lead to a discernible 
deficit in the SED (see ij |5.1[) . and the models presented 
to date are not compatible with large 40 - 50 AU gaps in 
massive disks which are still actively accreting. Further- 
more, due to shadowing, there are portions of the outer 
wall that will not be affected by the photoevaporative 
wind. 

Given the detection of over one hundred Jupiter- 
mass planets (Bu tler et al.l l2006h and about one dozen 
Neptu ne-mass planets (see references within lEndl et afl 
(2008)), it is natural to wonder if we can see evidence for 
the incipient stages of planet formation and particularly 
if pre-transitional disks around TTS constitute such evi- 
dence. At the moment, ALMA is the best suited facility 
to detect d isks around young planets within TTS disks 
(jWolf fc D 'Angclo 2005) and therefore search for planets 
within the gaps of pre-transitional disks. ALMA, along 
with JWST, will help us to better constrain the statis- 
tics of pre-transitional disks and their evolution. High- 
resolution near-IR intcrfcrometry can also examine the 
inner disk and reveal the spatial distribution of the opti- 
cally thin dust within the disk gap. 



6. SUMMARY & CONCLUSIONS 

In this paper we presented spectral type measurements 
as well as veiling and near-infrared excess measurements 
for the diskless LkCa 14, the pre-transitional disks of 
LkCa 15, UX Tau A, and Rox 44, and the transitional 
disks of DM Tau and GM Aur. Using near-IR SpeX 



spectra from 1-5 /im, we extracted the near-IR excess 
continua of our pre-transitional disk sample and fit their 
excess emission with optically thick inner disk walls, 
supporting the interpretation that these objects contain 
gaps in their disks. We modeled the broad-band SEDs 
of these pre-transitional disks and measure gap sizes of 
~40-70 AU. In the case of LkCa 15 and Rox 44 there is 
~10~ n M Q of ISM-sized optically thin dust within the 
gap. 

The near-infrared emission of our pre-transitional disks 
differs significantly from that of the transitional disks of 
DM Tau and GM Aur. DM Tau has no veiling and no 
excess emission in the near-IR, indicating that its disk 
hole is relatively empty of small dust grains. GM Aur 
has a small amount of veiling and displays a fiat excess 
above the stellar photosphere in the near-IR, which can 
be fit by a model of optically thin dust emission. This is 
in contrast to the pre-transitional disks which have large, 
blackbody-like near-IR excess continua that can be fit 
with models of a wall located at the dust sublimation 
radius. 

We also studied the effects of shadowing of the outer 
disk wall by the inner disk. We found that when the 
finite size of the star is taken into account, a significant 
portion of the outer wall is either in the prenumbra or 
completely out of the shadow of the inner wall. The 
predicted height of the wall is consistent with that of the 
outer disk, which is derived independently. 

Based on currently known disk clearing mechanisms, 
we propose that the gaps in pre-transitional disks are 
indicators of planet formation, making these disks a 
promising location for young planet searches. ALMA 
will be pivotal in extending the sample of known pre- 
transitional disks and has the potential to detect planets 
in these gaps. Near-IR intcrfcrometry will play an impor- 
tant role in further understanding the innermost regions 
of these disks. Future studies of this class of objects may 
bring us a few steps closer to understanding the origin of 
our own solar system. 
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APPENDIX 
THE EFFECT OF DISK SHADOWING 
A point source does not cast a penumbra. It only casts an umbra with a size, measured from above the midplane, of 



Zpoint — %w 



(Al) 
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where z w is the height of the inner wall above the midplane, r w is the radial distance between the center of the star 
and the inner wall, and r is the radius of the outer wall as measured from the center of the star. 

However, since the star is a finite source, the size of the umbra will be much smaller. At the outer wall, the height 
over which the whole star can be seen and the height below which the star cannot be seen are calculated using the 
roots of a quadratic function. This function is constructed for the distance between the inner wall and the points 
where rays tangent to the upper and lower parts of the star and the upper part of the inner wall cross the horizontal 
plane (Figure [TJ]). Using dimensionless distances (with i?» as unity) 




f) + pb/[(l + f)(l-f)]> (A2) 



x- = {r w -* (1 + _)(!-_) + if)/[(l +)(!)], (A3) 

Zy = (r - r w )(z w /x~) + Zw, (A4) 

z~ = (r-r w )(z w /x~) + z w , (A5) 

and, consequently, the outer wall above z~ will be illuminated to some degree. 

Given that a significant portion of the outer wall is not fully illuminated by the star, we consider the effect of the 
shadow, both the umbra and penumbra, on the heating of the wall and hence the resulting emission. The solid angle 
subtended by the star as seen by any point P in the disk, with coordinates r and z, is 

fl*(r,z) = 2 \ siruj>d(j)d6 (A6) 

Jo Jo 

where <j> characterizes an annulus at the stellar surface, and 9 allows one to integrate over the annulus. The maximum 
value of 4> is 

4>max = SWl -1 (l/d) (A7) 

where 

d=(z 2 + r 2 ) 1 / 2 . (A8) 

The maximum value of 9 depends on the shadow. 

In Figure [15] we show the relevant angles used to define the integration limits of the solid angle of the star. The 
angle 7 is related to the height of the shadow. If <j) < 7 the whole annulus can be seen, and then 9 max = ir. However, 
if <fi > 7 then a fraction of the annulus is in the shadow. 

From the figure, we can see that 

6max = 7T-E (A9) 

and 

To find the angle (3 we use that 

and 

Then, j3 is given by 



h = sva/3cosE. (A10) 

sin (3 = c sin <j> (All) 



c = d cos 4> - y/l- (dsin^) 2 . (A12) 



sin (3 — d cos <j> sin <j> - [1 - (d sin </>) 2 ] 1/2 sin tj>. (A13) 

The other relevant angle 7, is given by 7 = A — B, and so 

7 = tan.- 1 (r/z) - tan -1 [(r -r w )/(z- z w )}. (A14) 

The height h is given by 

h = (d- cos /3) tan 7 (A15) 

and thus, 

cos-E = (d — cos/3) tan 7/ sin /3, (A16) 

which is used in Eq. lA9l to calculate 9 max , which in turn allows one to calculate the solid angle, f2», with Eq. lA6l 
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Table 1 

Log of SpeX Observations 



Target 


KJ 


Date 


LXD 


Exposure (s) 


SNR^ 


LkCa 14 
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2009 Feb 18 
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LkCa 15 


8.2 


2009 Dec 22 
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370 
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7.6 


2009 Dec 22 






2400 
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Rox 44 


7.6 


2009 Apr 07 






1800 


150 


GM Aur 


8.3 


2009 Dec 22 






2400 


300 


DM Tau 


9.5 


2009 Feb 16 






2400 


170 


1 Two-Micron All- 


-Sky Survey Point 


Source Catalog 


(2MASS; 


Skrutskie et a] 


.1120061) 











Table 2 

Characteristics of Sample 



Target 


Spectral 


*K 


T' 

wall 

(K) 1 


Disk 




Type 




Class 2 


LkCa 14 


K5.5 







diskless 


LkCa 15 


K3 


0.6 


1400 


pTD 


UX Tau A 


G8 


0.4 


1550 


pTD 


Rox 44 


K3 3 


0.8 


1200 


pTD 


GM Aur 


K5.5 


0.1 




TD 


DM Tau 


Ml. 5 







TD 



^waii re ^ ers to the temperature of the best-fit 
wall to the the excess ncar-IR emission, if present. 

Our targets are classified as diskless stars, stars 
with transitional disks (TD), and stars with pre- 
transitional disks (pTD). 

Value adopted from Bouvicr & Appciizcllcr 



Table 3 

Stellar and Model Properties 
of LkCa 15, UX Tau A, & Rox 44 



Stellar Properties 






LkCa 15 


UX Tau A 


Rox 44 


M» (M ) 




1.3 


1.5 


1.3 


R* (R©) 




1.6 


1.8 


1.6 


T. (K) 




4730 
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4730 


M (Mq yr" 


') 


3.3xl0~ 9 


l.lxlO" 8 


9.3xl0~ 9 


Inclination ( 


deg) 


42 1 
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Distance (pc) 
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A v 
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mer Wall 




a,max (/Jm) 
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TL„H (K) 
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4„« (AU) 




0.017 
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0.034 


K,n.U (AU) 




0.15 


0.15 


0.25 




Optically Thick h 


mer Disk 




^inner disk 


(AU) 


<0.19 


<0.21 


<0.4 


dinner disk 


(Mq) 


<2xl0" 4 


<6xl0~ 5 


<8xl0 -5 


Optically Thick Outer Wall 


&max (,/^ m j 




0.25 


0.25 


0.25 


TZall (K) 3 




95 


110 


120 


<all (AU) 




12.9 


13.8 


9.9 


R° „ (AU)- 


i 


58 


71 


36 


Optically Thick Outer Disk 


e 




0.001 


0.001 


0.01 


a 




0.0005 


0.004 


0.006 


M disk (m@; 


I 


0.1 


0.04 


0.03 


Simon ct al. 


(2000) 








z lAndrews et al.| d2009l) 







Since we take into account the varying illumination along 
the outer wall due to the shadowing, there will be a range of 
T° a! , . Here we report the temperature of the fully illuminated 
part of the wall and use this value to calculate R.£, ii. 
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Figure 1. FAST spectrum of DM Tau compared to spectral type standards. We measure a spectral type of M1.5±l for DM Tau. Some 
of the features used for spectral classification are marked. The strong H Q emission feature at ~6550 A is cut off for clarity. 
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Figure 2. FAST spectra of GM Aur, LkCa 14, and LkCa 15 compared to spectral type standards. We measure a spectral type of K5.5±l 
for both GM Aur and LkCa 14 and a spectral type of K3±l for LkCa 15. Some of the features used for spectral classification are marked. 
The strong H a emission features at ~6550 A in GM Aur and LkCa 15 have been trimmed. 
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Figure 3. FAST spectrum of UX Tau A compared to spectral type standards. We measure a spectral type of G8±2 for UX Tau A. Some 
of the features used for spectral classification are marked. 
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Figure 4. Near-infrared emission of the diskless T Tauri star LkCa 14 relative to a standard star. The top panel contains the SpeX 
spectrum of LkCa 14 (dark line) compared to a K5 standard (light line). LkCa 14 has absorption lines that are not veiled relative to the 
standard star and we derive a veiling value (r^) of for LkCa 14 between 2.20 - 2.28 /im. The bottom panel displays the near-infrared 
excess emission of LkCa 14. After we subtract the standard star from the target (bottom), we find that there is no excess emission above 
the photosphere of LkCa 14 (gray). The dotted line corresponds to no excess and the small discrepancy between ~ 3.5-4.2 ^m is likely 
due to spectral slope uncertainties. This indicates that SpeX spectra of dwarfs are good templates for the photospheres of young stars. 
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Figure 5. SpeX near-infrared emission of the pre-transitional disks of LkCa 15, UX Tau A, and Rox 44 and the transitional disk of 
GM Aur relative to a standard star. After scaling the targets by 1 + rji , we subtracted the standard star spectra (top, light line) from 
the veiling-scaled target spectra (top, dark line) to extract the shapes of the excess emission above the photosphere (bottom, gray). We 
measure an r^ of 0.6, 0.4, 0.8, and 0.1 for LkCa 15, UX Tau A, Rox 44, and GM Aur, respectively. Additional veiling measurements from 
ISalvk e t al. ( 200SJ at ~4.8 /im (open points) are shown. For the pre-transitional disks, we find that the near-infrared excess emission can 
be reproduced by the wall of an optically thick inner disk with a temperature of T^ ;i (solid line). In the case of GM Aur, we can reproduce 
the near-IR excess emission with an optically thin dust model (bottom, broken line) composed of sub-micron-sized dust located within 
the inner disk hole following Calvct et al. (2005). We note that GM Aur was too faint to extract any data at wavelengths greater than 
~4.2 (im and that, in each of the spectra, regions from ~1.2— 1.4 /im, ~1.8-1.9 /im, and ~2.8— 3.2 /im are strongly affected by uncertainties 
in the telluric subtraction. 
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Figure 6. Near-infrared emission of the transitional disk of DM Tau relative to a standard star. In the top panel, we compare the SpeX 
spectrum of DM Tau (solid line) to a Ml dwarf standard star (broken line). We derive rx=0 for DM Tau between 2.20 - 2.28 (im and find 
that after subtracting the standard from the target (bottom panel) there is also no excess at longer wavelengths. The dotted line in the 
bottom panel corresponds to zero excess. The small discrepancy between ~ 2.4—2.5 (im is likely due to spectral slope uncertainties. Note 
that DM Tau was too faint to extract any data at wavelengths greater than ~ 4.2 /im. 



22 



Espaillat et al. 



10 



-8 



U I I I I I 1 1 1 — I — I I I I III 



i i i 1 1 1 1 1 — i — i— rr 



T-q 



LkCa 15 




A(^m) 



Figure 7. Pre-transitional disk model of LkCa 15. The best-fit model to LkCa 15 (solid black line) consists of an inner optically thick 
disk which extends from the dust destruction radius to <0.26 AU and an outer disk from 58 - 300 AU. Within the inner 4 AU of the 
gap between the inner and outer disk, the re is ^lO -11 Ma o f ISM-sized optically thin dust. Separate model components are the stellar 
photosphere (magenta dotted line; Kcnyon & Hartmann 1995), the inner disk (gray short-long-dash), the outer disk (red dot-short-dash), 
and the optically thin small dust located within the disk gap (green long-dash). [See the electronic edition of the Journal for a color version 
of this figure.] 
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Figure 8. Solid angle of LkCa 15 (dotted line), UX Tau A (dashed line), and Rox 44 (dot- dashed line) as seen from different heights (z) 
at the radius of the outer wall. Q is given in units of Q*, the solid angle of the entire star. 
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Figure 9. Pre-transitional disk model of UX Tau A. UX Tau A's best-fitting model (solid black line) has an inner optically thick disk 
which extends from the dust destruction radius to <0.21 AU and an outer disk from 71 - 300 AU. The gap between the inner and outer 
disk is relatively empty of small dust grains. Separate model components are the same as used in Figure [7] [Sec the electronic edition of 
the Journal for a color version of this figure.] 
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Figure 10. Pre-transitional disk model of Rox 44. The best-fit model to Rox 44 (solid black line) consists of an inner optically thick 
disk which extends from the dust destruction radius to <0.4 AU and an outer disk from 36 - 300 AU. Within the inner 2 AU of the gap 
between the inner and outer disk, there is ~10 -11 Mq of ISM-sized optically thin dust. Separate model components are the same as used 
in Figure [7] [See the electronic edition of the Journal for a color version of this figure.] 
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Figure 11. Schematic illustration of the pre-transitional disks of LkCa 15, UX Tau A, and Rox 44 and the transitional disk of GM Aur. 
The black semicircles represent the central stars and the small, empty circles are ISM-sized optically thin dust. The vertical lines are the 
disk walls. For the wall of the outer disk, black represents the portion of the wall that is fully illuminated by the star, light gray corresponds 
to the part of the wall that is in the penumbra of the inner disk, and dark gray is the portion of the wall that is in the umbra. In each 
case, the surface of the outer disk is fully illuminated. The transitional disk of DM Tau, which is not shown here, has an inner disk hole 
similar to what is seen in GM Aur, but it is devoid of optically thin small dust. Structures here are not drawn to scale. 
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Figure 12. Simulation of a gapped disk compared to a typical full disk model. A pre-transitional disk with a gap extending from ~0.3 AU 
to 4 AU (solid line) has a SED that diverges from a full disk SED (dashed line) beyond 10 fira. This is due to the contribution of the 
illuminated portion of the outer wall in the gapped disk model. The dotted line corresponds to the stellar photosphere. 
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Figure 13. Comparison of the opacities adopted in this work to those from Pollack et al. (1994) and Draine & Lee (1984). 
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Figure 14. Limits of the shadow cast by a finite source. Above z£ the whole star is visible and below z„ the star is in the umbra. In the 
intermediate region, the surface of the star is partially visible (i.e. the wall is in the penumbra). z po i n t is the height of the umbra when 
the star is taken to be a point source. 
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Figure 15. Geometry of the star, inner wall, and outer wall used to define the integration limits of the solid angle subtended by the star. 



